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ABSTRACT 

We report results on multiband observations from radio to 7 -rays of the two radio-loud 
narrow-line Seyfert 1 (NLSyl) galaxies PKS 2004—447 and J1548-I-3511. Both sources 
show a core-jet structure on parsec scale, while they are unresolved at the arcsecond 
scale. The high core dominance and the high variability brightness temperature make 
these NLSyl galaxies good 7 -ray source candidates. Fermi-LAT detected 7 -ray emis¬ 
sion only from PKS 2004—447, with a 7 -ray luminosity comparable to that observed 
in blazars. No 7 -ray emission is observed for J1548-I-3511. Both sources are variable 
in X-rays. J1548-I-3511 shows a hardening of the spectrum during high activity states, 
while PKS 2004—447 has no spectral variability. A spectral steepening likely related 
to the soft excess is hinted below 2 keV for J1548-I-3511, while the X-ray spectra of 
PKS 2004—447 collected by XMM-Newton in 2012 are described by a single power-law 
without significant soft excess. No additional absorption above the Galactic column 
density or the presence of an Fe line is detected in the X-ray spectra of both sources. 

Key words: galaxies: active - gamma-rays: general - radio continuum: general - 
galaxies: Seyfert 


1 INTRODUCTION 

Narrow Line Seyfert 1 (NLSyl) galaxies represent a rare 
type of classical Seyfert galaxies. The strong featureless 
X-ray continuum and the strong high ionization lines shown 
by NLSyl are common in Seyfert 1. However, the optical 
permitted emission lines are narrow, i.e. more similar to 
Seyfert 2 galaxies, indicating a combination of properties 
from both types. Their optical spectra are characterized 
by narrow permitted lines (full width at half-maximum, 
FWHM $( 2000 km s“^), weak [Oiii]/A5007 emission line 

i fOiii]/H/3 < 3) and usu ally strong Fen emission lines 
Osterbrock fc Poggelll985l) . 

NLSyl are usually hosted in spiral galaxies, although some 
objects are associated with early-type SO galaxies (see e.g. 
Mrk705 and Mrk 1239; Markarian et al. 1989). 

Extreme characteristics are observed in the X-ray band, 
where strong and rapid variability is observed more fre¬ 
quently than in classical Seyfert 1. The X-ray spectrum of 
NLSyl is usually described by a power law which dominates 
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in the 2-10 keV energy range, and a soft X-ray excess 
at lower energies. The X-ray spectru m between 0.3 an d 
10 keV is steeper than in Seyfert 1 (ICrupe et al.l [ioi^ l. 
while t he photon index of the hard X-ray spectrum is 
similar llPanessa et ahlBoill ). The complex X-ray spectrum 
is interpreted in terms of either relativistic blurred disc 
reflecti on, or ionized/neut r al absorption cove ring the X-ray 
source (iFabian et al.l[2009l : iMiller et al.ll201^ . 

About 7 per cent of NLSyl are radio-loud (RL), with a 
smaller fraction (~2.5 per cent) exhibiting a high radio¬ 
loudness parameteiQ (R > 100; Komossa et al. 2006). In the 
radio band, RL-NLSyl usually show a compact morphology 
with a one-sided jet emerging from the bright core and 
extending up to a few parsecs. In some objects the radio 
emission extends on kpc scales (iRichards fc Lister! l2015l : 
iDoi et ai] I 2 OI 2 I : lAnton et al.l l2008h . The high values of 
both brightness temperature and core dominance suggest 


^ The radio-loudness is defined as the ratio bet ween the radio flux 
densi ty at 6 cm and the optical B-band flux iKellermann et al] 

Il989h . 
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the presence of non-thermal emission from relativistic jets 
llDoi et al.ll201lll . The measurement of superluminal motion 
in the RL-NLSyl SBS 0846+ 513 indicates Dopp ler boosting 
effects in relativistic jets (iD’Ammando et alJ l2013ah . It 
is worth noting that some evidence of pc-scale jet-like 
structure is also found in some radio-quiet (RQ) NLSyl, 
but the nature of the outflow is still under debate (e.g. 
iDoi et ^ |[2 w3). 

Strong evidence in favour of highly-relativistic jets 
in RL-NLSyl is the detection by the Large Area 
Telescope (LAT) on board the Fermi satellite of 7 - 
ray emission from 6 RL-NLSyl: PMN J 0948+0022, 
IH 03 23+342, PKS 1502 +036, PKS 2004-44 7 (lAbdo et all 
l2009l j. SBS 0846+ 513 llD’Ammando et ah l2012| j. and 
FBQS J1644+2619 llD’Ammando et alJ 1201^ 1. The 7 -ray 
emission is variable, showing flaring activity accompanied 
by a moderate hardening of the spectrum. The peculiar 
multiwavelength properties together with the 7 -ray flares 
make the RL-NLSyl more similar to blazars rather than 
classical Seyfert galaxies, at least at high energies. 
Relativistic jets produced by nuclear objects which are 
thought to be hosted mainly in spiral galaxies are some¬ 
how puzzling. RL-NLSyl are usually found at higher 
redshift than RQ-NLSyl and no optical morphological 
studies have been carried out so far, with the exception of 
IH 0323+342. The optical morphology of the host galaxy 
is co mpatible with either a one-armed spiral (IZhou et al.l 
l2007li or a ring-like structure produced by a recent merger 
i Leon-Tavares et al.l[20l3 : lAnt6n et al.l|2008l L 
In this paper we present results of a multiwavelength 
study, from radio to 7 -rays, of the RL-NLSyl J1548+3511 
and PKS 2004—447. J1548+3511 is a NLSyl at redshift 


a = 0.478 with a radio loudness R~110 a nd an esti¬ 
mate d black hole mass Mbh = IO^’^Mq llYuan et al.l 
l2008h . while PKS 2004—447, at redshift z = 0.24, has 
1710 < R < 6320 and the e stimated black hole mass is 
10 ®’^Mq (lOshlack et al.l bOOll L These sources have been 
selected on the basis of their high variability brightness 
temperature T’e.var = 10 ^® - lO^^K that is considered a 
good indicator for the presence of Doppler boosting in rela- 
tiv istic jets. Among t he RL-NLSyl presented in the sample 
by lYuan et al.l 1 I 2 OO 8 II . J1548+3511 is the only source with 
T’s.var as high as those found in 7 -ray-loud NLSyl. How¬ 
ever, no 7 -ray emission has been detected from this source 
so far. In this paper we aim at investigating differences 
and similarities between the 7 -ray emitter PKS 2004—447 
and the 7 -ray silent J1548+3511. The results for these 
two sources are then compared to what has been found 
for the other 7 -ray emitting RL-NLSyl in the literature, 
in order to investigate the peculiarity of this sub-class of 
objects. The information from the multiwavelength data of 
J1548+3511 and PKS 2004—447 is then used to model the 
spectral energy distribution (SED) of these two sources. 
The paper is organized as follows: in Sections 2, 3, 4 and 
5 we report the radio, 7 -ray Fermi-hAT, X-ray Swift and 
XMM-Newton data analysis. In Section 6 we present the 
results from the multiwavelength analysis. The discussion 
and the presentation of the SED modelling are given in 
Section 7, while a brief summary is in Section 8 . 
Throughout this paper, we assume the following cosmology: 
Flo = 71 km s“^ Mpc“^, Dm = 0.27 and Da = 0.73, in a 


flat Universe. 


2 RADIO DATA 
2.1 VLB A observations 

Multifrequency Very Long Baseline Array (VLBA) observa¬ 
tion (project code BO045) of J1548+3511 was carried out 
at 5, 8.4, and 15 GHz on 2013, January 2. The observation 
was performed in phase reference mode with a recording 
bandwidth of 16 MHz per channel in dual polarization 
at 512 Mbps data rate. The target source was observed 
for 45 min at 5 and 8.4 GHz, and for 75 min at 15 GHz, 
spread into several scans and cycling through frequencies 
and calibrators in order to improve the uv-cover&ge. The 
strong source 3C 345 was used as fringe finder and bandpass 
calibrator. The source J1602+3326 was used as phase 
calibrator. 

Data reduction was performed using the NRAO’s As¬ 
tronomical Image Processing System (AIPS). A priori 
amplitude calibration was derived using measurements 
of the system temperature and the antenna gains. The 
uncertainties on the amplitude calibration (ctc) were found 
to be approximately 5 per cent at 5 and 8.4 GHz, and about 
7 per cent at 15 GHz. Final images were produced after 
a number of phase self-calibration iterations (Fig. [T]). The 
1(7 noise (rms) level measured on the image plane is about 
0.08 mjy beam~^ at 5 and 8.4 GHz, and about 0.12 mjy 
beam“^ at 15 GHz. The restoring beam is 3.4x1.3 mas^, 
2.1x0 .8 mas^, and 1.1x0.4 mas^ at 5, 8.4 and 15 GHz, 
respectively. In addition to the full-resolution images, at 5 
and 8.4 GHz we produced “low-resolution” images using 
natural grid weighting and a maximum baseline of 100 MA. 
The low-resolution image at 8.4 GHz is presented in Fig. [21 
The restoring beam is 2.6x1.8 mas^. 

To study the parsec-scale structure of PKS 2004—447 we 
retrieved archival VLBA data at 1.4 GHz (project code 
BD050). The observation was performed on 1998, October 
13 with a recording bandwidth of 8 MHz per channel in 
dual polarization at 256 Mbps data rate. In this observing 
run PKS 2004—447 was observed as phase calibrator, for 
a total time of about 30 minutes. Seven VLBA antennas 
participated in the observing run. Due to the southern 
declination of the source the restoring beam is highly 
elongated in the North-South direction and is 16.0x4.1 
mas^. 

The data were reduced following the same procedure 
described for J1548+3511. The uncertainties on the ampli¬ 
tude calibration are CTc ~ 10 per cent. Final images were 
produced after a number of phase self-calibration iterations. 
Amplitude self-calibration was applied using a solution 
interval longer than the scan length to remove residual 
systematic errors at the end of the self-calibration process. 
The final image is presented in Fig. [3] The Icr noise level 
measured on the image plane is ~0.3 mJy beam“'^. 

The total flux density of each source was measured by 
using the AIPS verb TVSTAT, which performs an aperture 
integration over a selected region on the image plane. In 
case of bright and compact components, like the core, we 
used the task JMFIT, which performs a Gaussian ht to the 
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Figure 2. Low-resolution VLBA image of J1548+3511 at 8.4 
GHz. On the image we provide the peak flux density in mjy 
beam“^ and the first contour intensity (f.c.) in mJy beam~^, 
which corresponds to two times the off-source noise level mea¬ 
sured on the image plane. Contours increase by a factor of 2. The 
restoring beam is plotted in the bottom left corner. 


PKS 2004-447 1.4 GHz 



Figure 3. VLBA image at 1.4 GHz of PKS 2004—447. On the 
image we provide the peak flux density in mjy beam“^ and the 
first contour intensity (f.c.) in mjy beam~^, which corresponds 
to three times the off-source noise level measured on the image 
plane. Contours increase by a factor of 2. The restoring beam is 
plotted in the bottom left corner. 


source component on the image plane. For more extended 
sub-components, like jets, the flux density was measured 
by TVSTAT. The uncertainty in the flux density arises 
from both the calibration error CTc and the measurement 
error (t™. The latter represents the off-source rms noise 
level measured on the image plane and is related to the 
source size 0 obs normalized by the beam size 0 beam as 
(Jm = rmsx(0obs/0beam)^^^. The flux density error crs 
reported in Table [T] takes both uncertainties into account 
and corresponds to as = y/a'i + crj^. 


2.2 Archival VLA data 

To investigate possible flux density variability we retrieved 
Very Large Array (VLA) archival data for both J1548-I-3511 
and PKS 2004—447 (see Table [2j . Data were reduced fol¬ 
lowing the standard procedure implemented in the AIPS 
package. Primary calibrators are 3C 286 and 3C 48 and the 
uncertainties on the amplitude calibration are between 3 
per cent and 5 per cent. VLA images were obtained after a 
few phase-only self-calibration iterations. Both sources are 
unresolved on arcsecond scale. 

The errors on the VLA flux densities are dominated by the 
uncertainties on the amplitude calibration, being am ~ 0.1 
mJy beam”*^. 


3 FERMI-LAT DATA: SELECTION AND 
ANALYSIS 

The LAT on board the Fermi satellite is a 7 -ray telescope 
operating from 20 MeV to > 300 GeV, with a large peak ef¬ 
fective area (~ 8000cm^ for 1 GeV photons), an energy res- 
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Figure 4. LAT light curve of PKS 2004—447 in the 0.1-100 GeV 
energy range obtained from 2008 August 4 to 2014 August 4 
(MJD 54682—56873) with 3-month time bins. Downward arrows 
represent 2cr upper limits. 


olution typically ~10 per cent, and a field of view of about 
2.4 sr with single-photon angular resolution (68 per cent con¬ 
tainment radius) of 0? 6 at E = 1 GeV on-ax is. Details about 
the LAT are given by I Atwood et al] ll2009l) . 

The LAT data reported in this paper for PKS 2004—447 
and J1548-I-3511 were collected over the first 6 years of 
Fermi operation, from 2008 August 4 (MJD 54682) to 
2014 August 4 (MJD 56873). During this time, the Fermi 
observatory operated almost entirely in survey mode. The 
analysis was performed with the ScienceTools software 
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Figure 1. VLBA images of J1548+3511 at 5 GHz (left), 8.4 GHz (centre), and 15 GHz (right). On each image we provide the peak flux 
density in mjy beam“^ and the first contour intensity (f.c.) in mjy beam“^, which corresponds to three times the off-source noise level 
measured on the image plane. Contours increase by a factor of 2. The restoring beam is plotted in the bottom left corner of each panel. 


Table 1. VLBA flux density of J1548+3511 and PKS 2004—447. Column 1: source name; Column 2: source component; Columns 3, 4, 
5, and 6: flux density at 1.4, 5, 8.4, and 15 GHz, respectively. 


Freq. 

GHz 

Comp. 

5i,4 

mJy 

Ss 

mJy 

SsA 

mJy 

‘S'ls 

mJy 

J1548+3511 

c 

- 

12.9±0.6 

16.0±0.8 

20.8±1.5 


J 

- 

3.8±0.2 

0.8±0.1 

- 


JO 

- 

- 

1.5±0.1 

- 


J1 

- 

9.7±0.6 

5.1±0.3 

- 


Ext 

- 

11.0±0.8 

4.4±0.5 

- 


Tot 

- 

37.4±2.0 

27.8±1.4 

20.8itl.5 

PKS 2004-447 

C 

224ib22 

- 

- 

- 


J 

lOOiblO 

- 

- 

- 


J1 

210ib21 

- 

- 

- 


Tot 

534ib53 

- 

- 

- 


Table 2. Archival VLA data. Column 1: source name; Column 2: Observing frequency; Columns 3, 4: date and project code of the 
observation; Column 5: beam size; Column 6: observing time; Column 7: flux density. 


Source 

Freq 

GHz 

Date 

Code 

Beam 

arcsec^ 

Obs. time 

min 

Flux 

mJy 

J1548+3511 

1.4 

1999-04-13 

AL485 

49.31x46.17 

1 

136±4 


1.4 

2003-09-05 

AL595 

1.46x1.17 

1.2 

124±6 


4.8 

1994-05-07 

AK360 

0.74x0.39 

1.0 

74±2 


8.4 

1995-08-14 

AM484 

0.24x0.22 

0.6 

76±2 


8.4 

1999-04-13 

AL485 

8.46x7.65 

0.8 

54±2 

PKS 2004-447 

8.4 

1988-12-23 

APOOl 

1.29x0.19 

1 

330±17 


8.4 

1995-07-15 

AK394 

1.15x0.19 

4.5 

270±11 


8.4 

2005-03-16 

AK583 

4.63x0.56 

3 

219±11 


8.4 

2005-04-21 

AK583 

3.65x0.60 

6 

440±18 
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package version v9r33pC0. Only events belonging to the 
‘Source’ class were used. The time intervals when the rock¬ 
ing angle of the LAT was greater than 52° were rejected. 
In addition, a cut on the zenith angle (< 100°) was applied 
to reduce contamination from the Earth limb 7 rays, which 
are produced by cosmic rays interacting with the upper 
atmosphere. The spectral analysis was performed with the 
instrument response functions P7REP_S0URCE_V15 using an 
unbinned maximum-likelihood method implemented in the 
tool gtlike. Isotropic (‘iso_source_v05.txt’) and Galactic 
diffuse emission (‘glLiem_v05_revl.fit’) components were 
used to model the backgrounc0. The normalizations of both 
components were allowed to vary freely during the spectral 
htting. 

We analysed a region of interest of 10° radius centred at 
the location of our two targets. We evaluated the signifi¬ 
cance of the 7 -ray signal from the source by means of the 
maximum-likelihood test statistic TS = 2 x (logLi — logLo), 
where L is the likelihood of the data given the model with 
(Li) or with out (Lp) a po i nt so urce at the position of our 
target (e.g., iMattox et al.l Il996li . The source model used 
in gtlike includes all the p oint sources from the third 
Eermi-LAT catalogue (3FGL; lAcero et al.l l2015h that fall 
within 15° of the target. The spectra of these sources were 
parametrized by power-law, log-parabola, or exponential 
cut-off power-law model, as in the 3FGL catalogue. 

A hrst maximum-likelihood analysis was performed to 
remove from the model the sources having TS < 10 and/or 
a predicted number of counts based on the fitted model 
Apred < 1. A second maximum-likelihood analysis was 
performed on the updated source model. In the htting 
procedure, the normalization factors and the photon indices 
of the sources lying within 10 ° of the target were left as free 
parameters. For the sources located between 10° and 15° 
from the target, we kept the normalization and the photon 
index hxed to the values from the 3FGL catalogue. 

For PKS 2004—447, the ht with a power-law model to the 
data integrated over 72 months of Eermi-LAT operation 
in the 0.1-100 GeV energy range results in TS = 164, 
with an average hux of (1.59 ± 0.16) xl0“® ph cm~^ s“^, 
and a photon index F^, = 2.39 ± 0.06, corresponding to 
an energy hux of (1.09±0.11)x 10“^^ erg cm~^ s“^. Figure 
[ 4 ] shows the 7 -ray light curve of PKS 2004—447 for the 
period 2008 August 4-2014 August 4 using 3-month time 
bins. For each time bin, the photon index was frozen to the 
value resulting from the likelihood analysis over the whole 
period. The systematic uncertainty in the effective area 
llAckermann et al]|2012l l amounts to 10 per cent below 100 
MeV, decreasing linearly with the logarithm of energy to 
5 per cent between 316 MeV and 10 GeV, and increasing 
linearly with the logarithm of energy up to 15 per cent at 
1 TeV 0. Statistical errors dominate over the systematics 
uncertainty. All errors relative to 7 -ray data reported 
throughout the paper are statistical only. 

For J1548-I-3511, the Ht with a power-law model over 6 
years of Eermi-LAT operation results in TS = 1. The 
2(7 upper limit is 3.35x10“® ph cm“^ s“^ in the 0.1-100 


GeV energy range (assuming a photon index P-y = 2.4), 
corresponding to an energy flux <5.3x10“^^ erg cm“^ s“^. 


4 SWIFT DATA: OBSERVATIONS AND 

ANALYSIS 

The Swift satellite llGehrels et al.l[2004h performed twenty 
seven observations of PKS 2004—447 between 2011 May 15 
and 2014 March 16. J1548-I-3511 has not been observed by 
Swift so far. 

The observations of PKS 2004—447 were carried out with 
all three instrument s on board: the X-ray Telescope (XRT; 
iBurrows et al.ll200^. 0.2-10.0 key)i th e Ultraviolet/Optical 
Telescope ('UVOT: lRoming et al.ll2 005l. 170-600 nml and the 
Burst Alert Telescope (BAT: iBarthelmv et al.ll2005i . 15-150 
keV). The source was not present in the Swift BA T 70-month 
hard X-ray catalogue llBaumgartner et aklfeoisl l. 

The XRT data of PKS 2004—447 were processed with 
standard procedures (xrtpipeline vO.13.0), Altering, and 
screening criteria using the HEAsoft package (v6.15). The 
data were collected in photon counting mode for all the 
observations. The source count rate was low (< 0.5 counts 
s“^); thus pile-up correction was not required. Source events 
were extracted from a circular region with a radius of 20 
pixels (1 pixel = 2.36 arcsec), while background events were 
extracted from a circular region with radius of 50 pixels 
away from the source region and from other bright sources. 
Ancillary response files were generated with xrtmkarf, 
and account for different extraction regions, vignetting 
and point-spread function corrections. We used the spec¬ 
tral redistribution matrices in the Calibration data base 
maintained by HEASARC0. Short observations performed 
during the same month were summed in order to have 
enough statistics to obtain a good spectral fit. The spectra 
with low numbers of photons collected (< 200 counts) were 
rebinned with a min imum of 1 count per bin and the Cash 
statistic (ICashin~979ll was used. We fitted the spectra with 
an absorbed power-law using th e photoelectric absorption 
model tbabs (IWilms et al.ll 2000 ll . with a neutral hydrogen 
column density fixed to its Galactic value (3.17x 10^°cm“^; 
iKalberla et alJ[2005fl . The fit results are reported in Table [S] 


UVOT data of PKS2004—447 in the v, b, u, wl, m2, 
and w2 Alters were reduced with the task uvotsource 
included in the HEAsoft package v6.15 and the 20130118 
CALDB-UVOTA release. We extracted the source counts 
from a circle with 5 arcsec radius centred on the source 
and the background counts from a circle with 10 arcsec 
radius in a nearby empty region. The observed magnitudes 
are reported in Table |T] We converted the magnitudes into 
de-reddened flux densities by using the E(B-V) value of 
0.029 from Schlaflv fc Finkbeineil (1201 lli . the extinction 
laws by iGardelli et al.| ||l989|) and the magnitude-flux 
calibr ations by iBessell et al.l (Il998fl and iBreeveld et al.l 
( 120111 ) . 


^ http://fermi.gsfc.nasa.gov/ssc/data/analysis/software/ 

^ http://fermi.gsfc.nasa.gov/ssc/data/access/lat/Background\Models.html 

http://fermi.gsfc.nasa.gov/ssc/data/analysis/LAT_caveats.html ^ http: //heasarc .nasa.gov/ 
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Table 3. Log and fitting results of S'roi/t-XRT observations of PKS 2004—447 using a power-law model with an HI column density fixed 
to the Galactic value in the direction of the source. “Unabsorbed flux. 


Date 

MJD 

Date 

UT 

Net exposure time 

s 

Photon index 
Px 

Flux 0.3-10 keV“ 
10“^® erg cm“^ 

55686 

2011-05-05 

6778 

1.65 ± 0.25 

5.6 ±0.9 

55756/55761/55765 

2011-07-14/19/23 

7101 

1.62 ± 0.29 

5.5 ±0.7 

55783/55786 

2011-08-10/13 

2258 

1.15 ± 0.49 

6.4 ±2.1 

55809/55821 

2011-09-05/17 

8239 

1.38 ± 0.19 

11.8 ± 1.3 

55880 

2011-11-15 

7095 

2.02 ± 0.30 

4.4 ±0.7 

56000 

2012-03-14 

7306 

1.52 ± 0.30 

4.5 ±0.8 

56111/56120 

2012-07-3/12 

7120 

1.65 ± 0.24 

6.5 ±0.9 

56182/56192/56200 

2012-09-12/22/30 

17016 

1.63 ± 0.18 

5.9 ±0.6 

56480/56487 

2013-07-07/14 

22963 

1.42 ± 0.12 

10.4 ±0.7 

56562 

2013-09-27 

8361 

1.45 ± 0.20 

14.6 ± 1.5 

56578/56585/56592/56594 

2013-10-13/20/27/29 

12699 

1.41 ± 0.24 

8.7 ± 1.2 

56599/56618 

2013-11-03/19 

8241 

1.58 ± 0.18 

6.8 ±0.6 

56619 

2013-11-20 

12179 

1.43 ± 0.16 

15.5 ± 1.2 

56730 

2014-03-14 

7507 

1.68 ± 0.26 

6.7 ± 1.0 

56732 

2014-03-16 

9642 

1.48 ± 0.19 

9.0 ±1.0 


Results of the Swift-WOT observations of PKS 2004—447 

in magnitudes. 



Date (MJD) 

Date (UT) 

V 

b 

u 

wl 

m2 

w2 

55686 

2011-05-05 

18.48±0.18 

20.08±0.28 

19.17±0.19 

19.77±0.32 

> 19.86 

> 20.19 

55756 

2011-07-14 

> 18.98 

19.69±0.30 

19.31±0.28 

19.29±0.29 

> 19.40 

> 19.85 

55761 

2011-07-19 

18.85±0.35 

- 

> 18.89 

> 19.35 

> 19.19 

> 19.61 

55765 

2011-07-23 

> 19.33 

20.25±0.36 

19.25±0.20 

19.53±0.28 

> 19.70 

> 20.01 

55783 

2011-08-10 

- 

- 

- 

19.41±0.18 

- 

- 

55786 

2011-08-13 

18.48±0.28 

19.56±0.31 

19.29±0.33 

19.67±0.35 

> 19.87 

20.37±0.35 

55809 

2011-09-05 

>18.24 

>19.12 

18.59±0.32 

> 19.06 

> 18.98 

19.74±0.34 

55821 

2011-09-17 

18.28±0.25 

18.99±0.16 

18.77±0.18 

18.84±0.24 

> 19.30 

> 19.50 

55880 

2011-11-15 

> 18.61 

19.63±0.32 

> 19.25 

> 19.33 

> 19.40 

> 19.52 

56000 

2012-03-14 

> 18.94 

> 19.76 

> 19.41 

> 19.59 

> 19.70 

19.55±0.25 

56111 

2012-07-03 

- 

- 

- 

19.48±0.11 

- 

- 

56120 

2012-07-12 

- 

- 

- 

- 

19.94±0.35 

- 

56182 

2012-09-12 

18.45±0.24 

19.67±0.30 

18.91±0.22 

19.46±0.28 

19.62±0.14 

20.00±0.23 

56192 

2012-09-22 

> 18.90 

19.36±0.23 

18.72±0.19 

19.37±0.32 

19.70±0.27 

20.27±0.30 

56200 

2012-09-30 

- 

- 

19.10±0.09 

- 

- 

- 

56480 

2013-07-07 

- 

- 

18.87±0.09 

- 

- 

19.91±0.20 

56487 

2013-07-14 

- 

- 

18.56±0.07 

19.36±0.10 

- 

- 

56562 

2013-09-27 

- 

- 

- 

19.03±0.09 

19.28±0.11 

- 

56578 

2013-10-13 

- 

- 

- 

- 

19.81±0.23 

- 

56585 

2013-10-20 

- 

- 

- 

- 

- 

19.77±0.15 

56592 

2013-10-27 

- 

- 

18.69±0.21 

- 

- 

19.54±0.25 

56594 

2013-10-29 

- 

- 

- 

- 

19.91±0.23 

- 

56599 

2013-11-03 

- 

- 

- 

19.28±0.27 

- 

- 

56618 

2013-11-19 

- 

- 

- 

18.85±0.15 

- 

- 

56619 

2013-11-20 

- 

- 

18.30±0.06 

- 

- 

- 

56730 

2014-03-14 

- 

- 

- 

19.67±0.16 

> 19.56 

- 

56732 

2014-03-16 

- 

- 

18.59±0.11 

- 

- 

- 


5 XMM-NEWTON DATA: OBSERVATIONS AND 
ANALYSIS 

5.1 EPIC Observations and data reduction 

J1548-I-3511 was observed by XMM-Newton ll Jansen et al.l 
[20 o 3) on 2011 August 8 and 20 for a total observing time 
of 28 ks in both cases (observation IDs 0674320301 and 
0674320401). For this analysis we focus on the data from 
the EPIC pn, which was operated in the Full Window mode, 


with net exposure times of 19.4 ks and 22.9 ks for the two 
observations. 

XMM-Newton observed PKS 2004—447 on 2012 May 1 and 
October 18 for a total duration of 36 ks in both cases (obser¬ 
vation IDs 0694530101 and 0694530201). The EPIC pn and 
the EPIC MOS cameras (MOSl and MOS2) were operated 
in the full frame mode. The first observation has net expo¬ 
sure times of 18, 20, and 20 ks for the pn, MOSl and MOS2, 
respectively; the second observation has net exposure times 
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of 29, 35, and 35 ks. Differently from the first observation, 
in the second observation there is a significant increase of 
counts (~70 per cent) adding the data collected by MOSl 
and MOS2 to the pn data. Therefore, for the first observa¬ 
tion of PKS 2004—447 we focus on the data from the EPIC 
pn only, while for the second observation both the pn and 
MOS data were analysed. 

The data were reduced using the XMM-Newton Science 
Analysis System (SAS vl3.5.0), applying standard event se¬ 
lection and filtering. Neither observation was affected by 
background flaring. The source spectra were extracted from 
a circular region of radius 30 arcsec centred on the source, 
and the background from a nearby region on the same chip. 
To allow for fitting the spectra were binned to contain 
at least 20 counts per bin. All errors are quoted at the 90 
per cent confidence level for the parameter of interest (cor¬ 
responding to Ax^ = 2.7). 

5.2 X-ray spectral analysis 

5.2.1 J1548+3511 

The spectral fits were performed over the 0.3-10 keV energy 
range using XSPEC v.12.8.1. Although we present only the 
fits to the EPIC pn, the results were cross-checked for consis¬ 
tency with the MOS spectra. Galactic absorption of 2.19 x 
10^° cm~^ was included in all fits using the tbabs model. In 
the two weeks between the observations the unabsorbed flux 
increased from To,3-10 keV = (4.2±0.2) x 10“^®erg cm“^ s“^ 
to To, 3-10 kev = (5.3 ±0.2) X 10“^®erg cm“^ s“^. No strong 
variability was seen during the individual observations. 

The results of fitting a single power law and a broken power- 
law to the two spectra are summarised in Tableland shown 
in Fig. [5] While the single power-law model leaves positive 
residuals above 3 keV, the broken power-law is a good fit 
to the data. The improvement between the models is more 
significant in the second observation, when the source was 
brighter. In this model the spectral shape changes from a soft 
slope of Pi = 2.6±0.1 (2.5±0.1) below Ebreak = S.OIq ® keV 
(1.7T0.2 keV) to r 2 = 1.97°:! (r 2 = 1.7T0.2) above Ebreak 
for the first (second) observation. If we add to the models 
another neutral absorber at the redshift of the source the 
fits do not improve, showing that no intrinsic absorption is 
required. Furthermore, there is no detection of an Fe line, 
with 90 per cent upper limits on a narrow line at 6.4 keV of 
EW < 0.26 keV and EW < 0.13 keV for the first and second 
observation, respectively. Unfortunately, the data quality is 
not sufficient to obtain meaningful constraints on more com¬ 
plex models for the emission. 

5.2.2 PKS2004-447 

As for J1548±3511, the spectral fits were performed over the 
0.3—10 keV energy range using XSPEC v. 12.8.1. Although 
for the first observation we present only the fits to the EPIC 
pn, the results were cross-checked for consistency with the 
MOS spectra. Galactic absorption of 3.17 x 10^° cm“^ was 
included in all fits using the tbabs model. In the five months 
between the observations the unabsorbed flux increased 


from Eo. 3-10 keV = (4.5 ± 0.1) x 10“^® erg cm“^ s“^ to 
Eo. 3-10 kev = ( 6.8 ± 0.2) X 10“^® erg cm“^ s“^. No strong 
variability was seen during the individual observations. 

The results of fitting a single power law and a broken 
power-law to the two spectra are summarised in Table |6l 
No significant soft X-ray excess is observed below 2 keV. A 
simple power-law model (Fig.[6)l is sufficient to describe the 
data of the first observation. For the second observation 
the power-law is not an optimal fit (y^ = 273/241, Table 
but no significant improvement was obtained using a 
broken power-law model. We note a dip in the residuals of 
the second observation at ~0.7 keV. Adding to the model 
an absorption edge, the fit is slightly better (y^ = 263/239) 
with a threshold energy Ec = 0.777 q keV, compatible 
with the Ovii absorption, and a maximum optical depth 
T = 0.211q ^^. However, observations with better statistics 
are required to confirm this feature. 

In both observations, the photon index is Px ~1.7, 
consistent with a jet emission component. Similarly to 
J1548±3511, if we add to the models another neutral ab¬ 
sorber at the redshift of the source the fits do not improve, 
showing that no intrinsic absorption is required. Moreover, 
there is no detection of an Fe line, with 90 per cent upper 
limits on a narrow line at 6.4 keV of EW < 0.12 keV 
and EW < 0.05 keV for the first and second observation, 
respectively. 


5.3 Optical Monitor data 

The Optical Monitor (OM; iMason et al.l |2001^ on board 
XMM-Newton is a 30 cm telescope carrying six optical/UV 
filters, and two grisms. We used the SAS task omichain to 
reduce the data and the tasks omsource and omphotom to 
derive the source magnitude. 

The OM was operated during all the observations described 
in Section 5.1. J1548±3511 was observed twice in optical and 
UV bands. Average observed magnitudes for J1548±3511 
are reported in Table [T] No significant change of activity 
was observed for J1548±3511 between the two OM observa¬ 
tions. PKS 2004—447 was observed by OM in u band, with 
observed magnitudes u = 19.20 ± 0.03 and u = 18.92 ± 0.03 
for 2012 May 1 and October 18, respectively. Therefore, on 
a half-year time-scale a difference of ~0.3 mag was observed 
in optical for PKS 2004-447. 

6 RESULTS 

6.1 The NLSyl J1548±3511 

6.1.1 Radio properties 

In the radio band, the NLSyl J1548±3511 shows a core-jet 
structure with a total angular size of about 70 mas corre¬ 
sponding to a linear size of 420 pc (Fig. [!}. Component C 
accounts for the majority of the VLBA flux density, from 
about 34 per cent at 5 GHz up to 100 per cent at 15 GHz. 
Its inverted spectrum {a ~ —0.4; Sv oc indicates that is 
the source core. From the core region a one-sided jet emerges 
with a position angle of about 10° and bends to the east (po¬ 
sition angle of about 30°) at a distance of ~40 mas (240 pc) 


http://xmm.esac.esa.int/external/xnim_user_support/documentatlciahsafehesgdBSp/f»±ere the component J1 is observed. The 
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Energy (keV) Energy (keV) 

Figure 5. EPIC pn spectra and residuals of J1548+3511 fitted with a power law (upper panel) and a broken power-law (lower panel). 
The first observation is shown to the left and the second one to the right. 


Table 5. Fits to the 0.3-10 keV EPIC pn spectra of J1548+3511. Galactic absorption was included in all fits. 


Model 

Parameter 

Value Obs 1 

Value Obs 2 

Power law 

r 

2.47 ±0.06 

2.31 ±0.05 


Norm 

7.8 ±0.2 X 10-5 

9.0 ±0.2 X 10-5 


X^/d.o.f. 

129/129 

210/164 

Broken power-law 

Pi 

9 ^7+0-09 

^•O<_o.o8 

^•^y_0.08 


^break (keV) 

2 q+0.8 

1 -T-t-O.S 
'-0.3 


r2 

1 q+0-3 
-*-•^-0.4 

1.7±0.2 


Norm 

7.5t|^-^ X 10-5 

8.4lo4 X 10-5 


X^/d.o.f. 

110/127 

154/162 


higher resolution provided by 8.4-GHz data allows us to re¬ 
solve the innermost part of the jet into two compact compo¬ 
nents (labelled J and JO in Fig.[T}, which are likely jet knots. 
An extended low-surface brightness structure (labelled Ext 
in Figs. [T] and [2]) is observed at ~60 mas (360 pc) from the 
core. At 15 GHz only the core component is detected. 

The source is unresolved on the arcsecond scale sampled 
by VLA images. The similar flux density (~141 mjy) re¬ 
ported in the NVSS J Condon et ahlUoOlii ) and the FIRST 


llBecker et al.l[l995h indicates that no extended emission on 
arcsecond scale is present. On the other hand, VLBA obser¬ 
vations at 5 and 8.4 GHz could recover only about 50 per 
cent or less of the VLA flux density. This may be due to a 
combination of both variability and the presence of extended 
jet structure that cannot be imaged by the short baselines 
of the VLBA. The spectral index computed using VLA data 
between 1.4 and 8.4 GHz results in a moderately flat spec¬ 
trum with a ~ 0.3-0.4. However, the spectral index values 
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Table 6. Fits to the 0.3—10 keV EPIC spectra of PKS 2004—447. Galactic absorption was included in all fits. 


Model 

Parameter 

Value Obs 1 

Value Obs 2 

Power law 

r 

1.72 ±0.05 

1.69 ±0.03 


Norm 

7.0 ±0.2 X 10“® 

10.0 ±0.3 X 10-® 


X^/d.o.f. 

97/102 

273/241 

Broken power-law 

Fi 

1.76 ±0.07 

1.71 ±0.03 


^break (keV) 

2 Q-h2.0 

q q-|-1.2 
'^•'^-0.6 


r2 

1.63 ±0.19 

1.51 ±0.14 


Norm 

7.0 X 10“® 

10.4 ±0.3 X 10“® 


X^/d.o.f. 

96/100 

268/239 


Table 7. Results of the XMM-OM observations of J1548+3511 in magnitudes. 


Date (MJD) 

Date (UT) 

V 

b 

u 

wl 

m2 

w2 

55781 

55793 

2011-08-08 

2011-08-20 

18.10±0.10 

17.99±0.07 

18.30±0.03 

17.37±0.02 

17.28±0.02 

16.99±0.02 

16.99±0.02 

16.62±0.08 

16.68±0.08 

16.77±0.07 

16.78±0.09 


are strongly subject to the flux density variability observed 
in this source. 

A small part of the flux density from the extended low- 
surface brightness structure is recovered in the low resolution 
image at 8.4 (Fig.[2j, where the total flux density is ~28 mjy 
instead of 21 mJy measured on the high-resolution image. 
No significant difference between the flux density measured 
in low- and high-resolution images is found at 15 GHz. 

We investigated the source variability by the analysis of 
archival VLA data. Although the data sets considered are 
not homogeneous (i.e. different VLA configurations), the 
lack of extended emission on arcsecond scale implies that 
variation in the flux density should be intrinsic to the source 
and not related to the lack of short baselines in the extended 
VLA configurations. Evidence for intrinsic variability comes 
from observations at 8.4 GHz, where the highest flux den¬ 
sity was measured when the VLA was in the most extended 
configuration (Table [2]). 

For ea ch frequency we comp uted the variability index V fol¬ 
lowing [H^^ttZ^eL^I] ||200^ : 


V = 


(5'max Umax) (‘S'min + <7min ) 
('S'max f^max ) + (S min -b ^min ) 


( 1 ) 


where Smax and Smin are the maximum and minimum flux 
density, whereas Umax and Umin are their associated errors, 
respectively. 

In addition to the VLA flux density reported in Table [21 
at 1.4 GHz we considered the values from the FIRST and 
NVSS, while at 5 GHz we considered the values reported 
in the 87GB catalogue ([Gregory fc Cond ^ 1991^ and i n 
the second MIT-Green Bank survey ( Langston et aI.l|l99Ct l. 
VLBA flux densities were not taken into account due to 
the possible missing flux from extended jet structures on 
parsec scales. From Eq. [T] we found that V is 1 per cent, 
11 per cent and 14 per cent, at 1.4, 5, and 8.4 GHz respec¬ 
tively, indicating larger variability at higher frequencies as 
usually found in blazars. The low variability {V = 0.01) 
estimated at 1.4 GHz is comparable to the uncertainties. 


We note that the flux density variability may be under¬ 
estimated due to the poor time sampling of the observations. 


6.1.8 X-ray properties 

In the two weeks between the observations the source 
brightened by about 25 per cent. The XMM-Newton spec¬ 
tra of J1548-1-3511 are well fitted by a broken power-law 
model, with a possible Seyfert component below ~2 keV 
and a jet component dominating at higher energies. The 
low-energy component, with a steep photon index of ~2.5, 
may be associated with the soft X-ray excess. On the other 
hand, the relatively hard X-ray spectrum above the energy 
break (r 2 = 1.7-1.9) may suggest a significant contribution 
of inverse Compton radiation from a relativistic jet. 


6.1.3 Optical and UV properties 

No significant change of activity was observed for 
J1548-I-3511 in the optical and UV bands between the two 
AMM-OM observations performed two weeks apart (Table 

0 . 

6.2 The NLSyl PKS 2004-447 

6.2.1 Radio properties 

The radio source PKS 2004—447 has a core-jet structure 
with a total angular size of about 40 mas, which corresponds 
to ~150 pc at the redshift of the source. The radio emission 
is dominated by the source core, labelled C in Fig. [31 which 
accounts for ~42 per cent of the total flux density at 1.4 
GHz. The jet structure emerges from the core component 
with a position angle of about —90°, then at 20 mas (~75 
pc) it bends to a position angle of about —60°. The jet 
structure is resolved into two subcomponents, J and Jl, 
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Figure 6. EPIC spectra and residuals of PKS 2004—447 fitted 
with a power law. For the first observation (upper panel) only 
EPIC pn data were considered, while for the second (bottom 
panel) EPIC pn, MOSl, and MOS2 were considered in the fit. 


which are enshrouded by diffuse emission. The lack of 
multifrequency observations does not allow us to study the 
spectral index distribution across the source. 

The source is unresolved in VLA images at 8.4 GHz, 
in agreemen t with previous studies at arcsecond-scale 
resolution bv iGallo et al.l (l2006l) . 

Archival VLA observations at 8.4 GHz point out some 
level of flux density variability (Table [2|). From Eq. [T] we 
computed the variability index for PKS 2004—447, which 
turns out to be 27 per cent, consistent with the flux density 
varia bility derived from the Ceduna observations at 6.65 
GHz dGallo et ahllioo^) . 


6.2.2 'y-ray propeHies 

During the first six years of Fermi-LAT observations, the 
0.1-100 GeV averaged flux is ~1.6xl0“® ph cm“^ s“^. The 
LAT light curve indicates variable 7 -ray emission with flux 
ranging between (1.3-4.2)x 10“® ph cm“^ s“^, interleaved 
by periods of low activity, when the source is not detected 


by Fermi-LAT (Fig. |4]). No 7 -ray flares from this source 
have been detected so far. 


6.2.3 X-ray properties 

The X-ray light curve collected by Swift-XRT indicates 
significant flux variability between 2011 and 2014, ranging 
between (5-16) x 10 “^^^ erg cm“^ s“^ (Fig- O- In 2011 
September the increase of the X-ray flux occurs when the 
7 -ray emission is in a maximum, suggesting a possible 
correlation between the emission in these two energy 
bands. On the other hand, between 2013 September and 
November the X-ray light curve shows two high-activity 
episodes in which the flux is F 0 . 3 —lokev ~ 15x10 erg 
cm“^ s“^, interleaved by low-activity state. During the 
X-ray high-activity state in 2013 November, the 7 -ray flux 
is roughly 2.6 times the average value, whereas in 2013 
August-October the source is not detected in 7 -rays. No 
significant X-ray photon index variability is observed (Fig. 
0 . 

In the five months between the XMM-Newton observations 
the source brightened by about 35 per cent. The photon 
index derived for the two XMM-Newton observations is in 
good agreement with the Swift-XMT results, an d is slightly 
softer than the value obtained in 2001 by iGallo et al.l 
(l2006l) . This may be due to the higher flux observed in 
2001 (Tc.B-iofeeV = 1-5 X 10“^^ erg cm“^ s“^) than that 
observed in the 2012 XMM-Newton observations. In the 
2012 XMM-Newton observations, there is no evidence of 
soft X-ray excess. 


6 . 2.4 Optical and UV properties 

On monthly time-scales a difference of ~0.3 mag was ob¬ 
served for PKS 2004—447 in u band by XMM-OM.. Dur¬ 
ing the Swift-WOTt observations the change of magnitudes 
spanned in the various band is about 0.6, 1.3, 1.0, 0.9, 0.7, 
0.8 going from the v to w2 filter (Table[4)). A first peak of the 
optical/UV activity was observed on 2011 September 5, dur¬ 
ing a high X-ray and 7 -ray activity period. A second peak of 
activity was observed in the u and wl bands on 2013 Novem¬ 
ber 19 and 20. At that time, the maximum X-ray flux was 
observed together with a high 7 -ray flux level. A contem¬ 
poraneous increase of activity between the optical-UV and 
the X-ray and 7 -ray bands indicates that the jet emission 
is dominant also in optical and UV, in agreem ent with the 
lack o f a significant disc emission suggested bv lAbdo et al.l 
(l2009l) . 

7 DISCUSSION 

7.1 Relativistic jets in RL-NLSyl 

High-energy emission has been detected in RL-NLSyl. 
On the other hand, no 7 -ray emission has been found in 
RQ-NLSyl, suggesting a possible intrinsically different 
nature between these two sub-populations. Different 7 -ray 
properties have been observed in the six RL-NLSyl de¬ 
tected by Fermi-LAT. Three objects, PMN J0948+0022, 
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Figure 7. Multifrequency light curve of PKS 2004—447 for the period between 2011 May and 2014 August. The LAT q-ray light curve 
(in units of 10~® ph cm~^ s“^ for 0.1 < E < 100 GeV; upper panel), the Swift X-ray light curve (in units of 10“^® erg cm~^ s“^ for 
0.3 < E < 10 keV; second panel from top), the X-ray photon index trend {third panel), the Swift-WOTt light curve in the optical and 
UV filters {v, b, u, ml, m2, w2, in units of mjy; fourth to ninth panel). Arrows refer to 2cr and Scr upper limits on the source flux for 
LAT and Swift-WOT measurements, respectively. 


SBS 0846+513, and IH 0323+342 show 7 -ray flares, reach¬ 
ing a n intrinsic apparen t lumin o sity a s high as 10 "^® erg 


(iD’Ammando et al.l l2015bl . 2012 f). i.e. comparabl e 


___ n parabl e 

to those shown by FSRQ (e.g. lAckermann et alT 1201 j l. 
On the other hand, PKS 2004—447, PKS 1502+036, and 
FBQS J1644+2619 have not shown strong flares so far. 
During the first 6 years of Fermi operation, the apparent 
luminosity of PKS 2004—447 ranges between L~f = (1.3- 
4.2) X10“*® erg s“^. A similar behaviour was found for 
PKS 1502-1 -036 with a slightly larger l uminosity {Lj ~ 10“^® 
erg s“^, ID’Ammando etlil] l2013bl l. FBQS J1644+2619 
had an intermediate behaviour with hi gh activity episodes 
interle aved by long quiescent periods (ID’Ammando et al.l 
l2015all . 

The RL-NLSyl J1548+3511 is not detected in 7 -rays and 
the upper limit to the apparent isotropic luminosity is 


Lj < 1.7 X lO'^® erg s“^. 

The detection of 7 -ray emission in a handful of RL-NLSyl 
proves the presence of relativistic jets in this peculiar 
sub-class of active galactic nuclei (AGN). In the photon 
index versus 7 -ray luminosity plane, the 7 -ray loud NLSyl 
are located In the low-luminosity tail of FSRQ distribution, 
where also the 7 -ray emitting steep-spectrum radio q uasars 
are found (lAckermann et alll2015l : lAbdo et aI.ll2010a L 
A jet component contribution is likely observed above 2 
keV in the RL-NLSyl J1548+3511, where the X-ray photon 
index is Fx ~ 1.7-1.9. Below 2 keV the spectrum is softer, 
compatible with the presence of soft X -ray excess usually 
observed in NLSyl (iGrupe et al.l l20ld') , as well as in the 
7 -ray NLSyl PM N J0948+0022 J^mmando et ID l2014ll 
and IH 0323+342 (IPaliva et al.l |2014|L An indi cation of a 
weak soft excess was reported by Gallo et al.l (I 2 OO 6 II for 
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PKS 2004-447 in 2004. The X-ray spectra of PKS 2004-447 
observed in 2012 are well reproduced by a single power law 
with hard photon index Px ~ 1.7 and no significant soft 
X-ray excess is needed below 2 keV. This may be due to 
the relatively low statistics of the 2012 observations. The 
study of the multiwavelength variability of PKS 2004—447 
indicates that during the high activity state observed 
in 7 -rays between 2011 August and October, when the 
7 -ray flux was a factor of 2.6 times the average value, also 
the X-ray, UV and optical emission reached a maximum, 
suggesting a common origin for the multiband variability. 
It is worth mentioning that not all the episodes of flux 
increase in X-rays/UV are associated with a high activity 
state in 7 -rays (Fig. [7]), like in the case of the high X-ray 
activity observed in 2013 September. Such X-ray/UV flares 
with no obvious counterpart in other bands were observed 
in the RL-NLSyl PMN J0948-I-0022 (e.g. |FTOclnrn_^£_al 
201^ as well as in many blazars (e.g., 3C 279. lAbdo et~ 


No systematic studies of complete samples of NLSyl 
have been carried out at Very High Energy (VHE) so 
far. PKS 2004—447 was observe d by the High Energ y 
Stereoscopic System (H.E.S.S., lAharonian et al] I2OO6II . 
No detection was obtained and the esti mated upper limit 
was 0.9 per cent of the Crab Units llAbramowski et aP 
l2014h . Upper limits of 10 per cent and 1.9 per cent of 
the Crab Uni ts at VHE were o btained for the NLSyl 
IH 0323-f342 dFalcone et al.l l2004ll and PMN J0948-f0022 
llD’Ammando et^iL 2015lJl . respectively. The lack of VHE 
detection may be due either to an intrinsically soft 7 -ray 
spectrum, or to 7-7 absorption within the source, or to 
pair production from 7 -ray photons of the source and the 
infrared photons from the extragalactic background light 
(EBL), although the latter scenario is disfavoured by the 
fact that the redshift o f the most distant FSRQ detecte d 
at VHE, PKS 1441-1-25 (iMirzovan ll2015l : lMukheriee Il2015ll . 
is higher {z = 0.939) than the redshifts of the RL-NLSyl 
investigated so far. 


T.2 Physical properties 

Milliarcsecond resolution observations are a fundamental re¬ 
quirement for describing the morphology and understanding 
the physical properties of RL-NLSyl. Given their compact¬ 
ness on arcsecond scale, the high angular resolution provided 
by VLBA observations allows us to investigate the presence 
of a jet structure emerging from the core region and to con¬ 
strain the physical characteristics of the core emission with¬ 
out substantial contamination from the jet. 

We computed the brightness temperature Tb of the core 
component by using: 


Tb 


1 g(^) (cy 

2k Q \v) 


( 2 ) 


where k is the Boltzmann constant, fl is the solid angle of 
the emitting regions, S{i') is the source-frame flux density 
at the observed frequency, and c is the speed of light. We 
computed the source-frame flux density as S{u) = Sohs{v) x 
(1 + zA~°‘, where Sobsii^) is the observer-frame flux density, 
z is the source redshift, and a is the radio spectral index 


that we assume to be equal to 0 for the self-absorbed core 
component. The solid angle is given by: 

H = — ^min^maj (3) 

where 0 maj and 0 min are the major axis and minor axis, 
respectively. In both sources the core region is unresolved 
by our VLBA observations and the major and minor axes 
must be considered as upper limits. This ends up in a lower 
limit to the brightness temperature. 

If in Eq. [2] we consider the values derived from the 
VLBA image at 15 GHz, we obtain a brightness temper¬ 
ature Tb > 4.4 X 10® K for J1548-I-3511. In the case of 
PKS 2004—447 we computed the core brightness temper¬ 
ature making use of the values derived from the VLBA 
image at 1.4 GHz, and we obtained Tb > 2 x 10^® K. 


For PKS 2004—447, the availability of two VLA ob¬ 
servations separated by about one month allowed us to 
estimate the rest-frame variability brightness tempera- 
ture, Tb var by using the f lux density variability. Following 
iD’Amrnando etlil] (l2013bll we computed T^ by: 

T' 2 \^S\Dl 

nk + z)^+‘^ 


where |AS| is the flux density variation. At is the time lag 
between the observations, and T>l is the luminosity distance. 
If in Eq. |4] we consider IASI = 221 mjy and At = 36 
days, i.e. the flux density variation measured between the 
last two VLA observations, we obtain T^ var ~ 1-7 x 10^‘^K. 
This high variability brightness tempera ture is similar to 
the value derived by ICallo et al.l (l2006h on the basis of 
Ceduna monitorin g campaign. In the case of J1548+3511, 
lYuan et al.l ll2008h estimated a variability brightness tem¬ 
perature Ta^var — 10 ^®K on the basis of the flux density 
variation measured in 207 days at 4.9 GHz. Assuming that 
such high values are due to Doppler boosting, we can esti¬ 
mate the variability Doppler factor (5var, by using: 


5 


var — 


TB.var \ 
Tint J 


(5) 


where Tnt is the intrinsic brightness temperature. Assuming 
a typi cal va lue Tnt = 5 x 10 ®®K, as derived by e.g. iReadheadI 
(Il994h and iHovatta et ^ ll2009l j. and a flat spectrum with 
a = 0, we obtain 5 = 15 and 5.8 for PKS 2004-447 
and J1548-I-3511, respectively. These values are in agree¬ 
ment with the range of variabili ty Doppler factor found fo r 
the RL-NLSyl SBS0846-f513 j D’Arnrn ando et ^12013^ 1. 
PKS 150 24-036 (iD’Am mand o et al.l 2013blj . as well as for 
blazars (ISavolainen et all 201Cllj . 

We estimated the ranges of viewing angles 8 and of the bulk 
velocity in terms of speed of light /3 from the jet/counter-jet 
brightness ratio. Assuming that the source has two symmet¬ 
rical jets of the same intrinsic power, we used the equation: 


_ / l-b/3cosg y+" 

Bcj V1 “ /3cos6 ') 


( 6 ) 


where Bj and B^j are the jet and counter-jet brightness, 
respectively. We prefer to compare the surface bright¬ 
ness instead of the flux density because the jet has a 
smooth structure without clear knots. The jet brightness 
for J1548-b3511 and PKS 2004-447 is 0.93 mJy beam"i 
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and 14.0 mjy beam“^, respectively. In the case of the 
counter-jet, which is not visible, we assumed an upper 
limit for the surface brightness that corresponds to 0.2 
mJy beam“^ and 0.6 mJy beam“^ for J1548-I-3511 and 
PKS 2004—447, respectively, i.e. la noise level measured 
on the image. From the brightness ratio estimated from 
Eq. [ 6 ] we obtain l3cos6 > 0.61 for J1548-I-3511 and 0.52 
for PKS 2004—447, implying that the minimum velocity is 
P > 0.6 and 0.5 and a maximum viewing angle 6 = 52° 
and 59°, respectively. These limits do not provide a tight 
constraint on the physical parameters of these objects. 


7.3 SED modelling 

We created an average SED for the two NLSyl studied 
here. The SED of PKS 2004—447 includes the 6 -year 
average Fermi-LAT spectrum, the XMM-Newton EPIC-pn 
data collected on 2012 October 18, and the Swift-UYOT 
data collected on 2012 September 12. In addition, we 
included the IR data collected by WISE on 2010 April 14 
and by Siding Spring Observatory in 2004 April 10 - 12, 
the radio VLA and VLBA data presen t ed in this paper, 
and the ATCA data from lOallo et al.l ll2006l j. The SED 
of J1548-I-3511 includes the upper limit estimated over 6 
years of Fermi observations, the XMM-Newton (EPIC-pn 
and OM) data collected on 2011 August 8 , the IR data 
collected by WISE on 2010 January 30 and by 2MASS on 
1998 April 3, and the radio VLA data presented here. The 
multiwavelength data are not simultaneous. 

We modelled the two SED with a combination of syn¬ 
chrotron, synchrotron self-Compton (SSC), and external 
Compton (EC) non-thermal emission. The synchrotron 
component considered is self-absorbed below 10^^ Hz and 
thus cannot reproduce the radio emission. This emission 
is likely from th e superpositi on of multiple self-absorbed 
jet components JKonigll Il981^ . We also included thermal 
emission by an accreti on disc and dust to rus. The modelling 
detail s can be found in iFinke et al.l l|2008l 'l and lPermer et aP 
l|2009h . Additionally, a soft excess was observed in the X-ray 
spectrum of J1548-I-3511. We note that the origin of this 
soft X-ray emission is uncertain. In order to account for 
this feature, as a possible origin, we included emission from 
the disc, which is Compton scattered by an optically thin 
thermal plasma near the accretion disc (i.e., a corona). Thi s 
was done using the routine “SIMPL” (jSteiner et ahllMOSh . 
This routine has two free parameters: the fraction of disc 
photons scattered by the corona (/sc), and the power-law 
photon index of the scattered coronal emission (Fsc). The 
mass of th e BH was chosen to be the same as the one 
reported in iFoschini et ahl l|2015h . 

The results of the modelling are reported in Table [ 8 ] and 
Figure [ 8 ] (for a de scription of the model parameters see 
lDermeretal.1 [ 200 ^ 1 . This model assumes that the emitting 
region is outside the Broad Line Region, where dust torus 
photons are likely the seed photon source. This seed photon 
source was modelled as being an isotropic, monochromatic 
radiation source with dust parameters chosen to be consis¬ 
tent with the relation betw een inner radius , disc l uminosity, 
and dust temperature from iNenkova et al .1 (l2008l j. 

The Compton component of the SED of PKS 2004—447 
and J1548-I-3511 is modelled with an external Comp¬ 


ton scattering of dust toru s seed photons, as for the 


SED of SBS 0846-1-513 _ ( DWmrnundoeL^ _ 2013lJ l 


and PMN J0948-I-0022 l|D’Ammando et al.l l2015bl 'l. In 
PKS 2004—447 the IR data are not well fitted by the model. 
This may be due to the fact that the data considered in 
the SED are not simultaneous and some variability may 
be present (note th at these data are no t well-fitted by the 
modelling done by IPaliva et al.1 (120131 ') either). We note 
that the disc luminosity of PKS 2004—447 is particularly 
weak. This value is similar to the disc luminosity of 
SBS 0846-1-513, while it is significantly lower than that of 
PMN J0948+0022. The disc luminosity of PKS 2004-447 
is con sistent with the value estimated by iFoschini et al.1 
(l2015lj on the basis of the optical spectrum, and the lack 
of a blue bump. This is in co n trast to the modelling of 
PKS 2004—447 by IPaliva et al.l ll2013ll who used a much 
brighter dust component, which w as not consiste n t wit h 
the disc luminosity estimated by IFoschini et al.1 ll2015h . 
It is worth mentioning that the model parameters shown 
here are not unique, and other model parameters could 
reproduce the SED equally well. This is particularly true for 
J1548-I-3511, since without a LAT detection, its Compton 
component is not well-constrained. 


7.4 RL-NLSyl and the young radio source 
population 

The low black hole mass and the high accretion rate 
commonly estimated in NLSyl suggeste d that these objects 
may be in an early evolutionary stage ijCruDe et al.1 Il999l : 
lMathuill2000l '). Some RL-NLSyl have been proposed to be 
compact steep spectrum (CSS) radio sources, on the basis of 
their compactness and flat/inverted spectrum which turns 
over at a few hundred MHz or around one GHz (lYuan et al.1 
l2008h . CSS are powerful radio sources whose linear size is 
< 20 kpc. Due to their intrinsically compact size they are 
considered radio sources in an early evolutionary stage (see 
e.g. O’Dea 1998 for a review). Kinematic and radiative 
studies provided ages of 10 ^~^ years, strongly favouring the 
youth scenario (see e.g. IPolatidis &: Conwavl l2003l : iMurgial 
l2003h . 

Many CSS sources are hosted in galaxies that recently 
underwent major mergers, which may have triggered 
the onset of the radio emission . A s imilar scenario was 
suggested bv iLeon-Tavares et al.1 l|2014l ') for the RL-NLSyl 
IH 0323+342. A link between NLSy l and CSS was also 
proposed by ICaccianiga et all (l2014l ~l for the RL-NLSyl 
SDSS J143244.91+301435.3 on the basis of its compact 
size, absence of variability, and a steep radio spectrum. In 
this context the variable and 7 -ray emitting RL-NLSyl 
may be the aligned population of RL-NLSyl where the 
CSS properties are hidden by dominant boosting effects. 
However, at least three (IH 0323+342, PMN J0948+0022, 
and FBQS J1644+2619) out of the six RL-NLSyl detected 
in 7 -rays have extended structures with linear size of 
20-50 kpc (|Doi et al ]|2 qI 3), challenging this interpretation. 
Among the remaining three objects, PKS 2004—447 wa s 
suggested as a possible CSS source bv ICallo et al.l (l2006li . 
However, the X-ray spectra in CSS sources are typically 
highly obscured with c olumn density Nh > 10 ^^ cm“^ 
jTengstrand et al1l2009h . while no absorber in addition to 
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Table 8. Parameters used to model the SED. 


Parameter 

Symbol 

PKS 2004-447 

J1548+3511 

Redshift 

z 

0.24 

0.478 

Bulk Lorentz Factor 

r 

30 

30 

Doppler factor 


30 

30 

Magnetic Field [G] 

B 

0.5 

2.0 

Variability Timescale [s] 

t-i] 

10® 

10® 

Comoving radius of blob [cm] 

K 

7.3xl0i® 

6.1 X lO^® 

Low-Energy Electron Spectral Index 

Pi 

2.5 

2.5 

High-Energy Electron Spectral Index 

p2 

3.8 

3.6 

Minimum Electron Lorentz Factor 

l' ■ 

' min 

1.0 

3.0 

Break Electron Lorentz Factor 

'T'brk 

6.0 X 10^ 

1.0 X 10^ 

Maximum Electron Lorentz Factor 

^max 

4.0 X 103 

1.0 X 10^ 

Black hole mass [Mq] 

A^bh 

4.3 X 10® 

8.3 X 10'^ 

Disc luminosity [erg s“^] 

-^disc 

1.8 X lO'^^ 

1.4 X 10"^® 

Inner disc radius [i?g] 

Rin 

6.0 

2.0 

Outer disc radius [Rg] 

Rout 

200 

200 

Fraction of disc photons scattered by corona 

/sc 

0 

0.08 

Corona spectral index 

Psc 

N/A 

2.6 

Seed photon source energy density [erg cm~^] 

^seed 

8.3 X 10“® 

2.7 X 10“® 

Seed photon source photon energy 

^seed 

5 X 10“’^ 

5 X lO-’’ 

Dust torus luminosity [erg s“^] 

-^dust 

7.5 X 10^® 

1.9 X 10** 

Dust torus radius [cm] 

-Rdust 

6.5 X lOi® 

4.1 X lOi® 

Dust temperature [K] 

Tdust 

1000 

1000 

Jet Power in Magnetic Field [erg s~^] 


8.9 X 10"*® 

1.0 X 10'^’^ 

Jet Power in Electrons [erg s~^] 

Pi,e 

2.7 X 10-*"* 

6.6 X 10"^® 




Figure 8. SED data (circles) and model fit (solid curve) of J1548H-3511 (left) and PKS 2004—447 (right). The model components are 
shown as dashed curves. Multiwavelength data are not simultaneous. 
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the Galactic one is needed for modelling the X-ray spectrum 
of PKS 2004-447. 

The information available so far is not enough to firmly 
link the RL-NLSyl to the young radio source popu¬ 
lation. Statistical multiwavelength studies on a large 
sample of RL-NLSyl are required for investigating a pos¬ 
sible connection between these sub-class of radio-loud AGN. 


8 CONCLUSIONS 

We presented results on a multiwavelength study, from radio 
to 7 -rays, of the RL-NLSyl J1548-I-3511 and PKS 2004—447. 
The conclusions from this investigation can be summarized 
as follows: 

• PKS 2004—447 is detected in 7 -rays by LAT with 
an average flux between 0.1 and 100 GeV of ~1.6xl0~® 
phcm“^ s“^, corresponding to a luminosity of ^ 1 . 6 x 10 “*® 
erg s“^, which is comparable to the values found in the 
other 7 -ray emitting NLSyl. No strong flares have been 
detected so far. 

• J1548-I-3511 has not been detected in 7 -rays by LAT 
during the first 6 years of observations. The upper limit to 
the luminosity is 1.7x10“*® erg/s. 

• Both sources have a clear core-jet structure on parsec 
scales. The majority of the radio emission comes from the 
core component. On arcsecond scale the radio structure is 
unresolved. 

• PKS 2004—447 shows significant variability from radio 
to 7 -rays. In particular, at the end of 2011 and at the end of 
2013 the high activity state observed in 7 -rays is simultane¬ 
ous to a local maximum in X-rays, UV and optical bands, 
suggesting a common origin. The X-ray spectra collected 
by XMM-Newton in 2012 from 0.3 to 10 keV are fitted 
by a single power law. No significant X-ray photon index 
variability is observed in the period considered in this paper. 

• The X-ray spectrum of J1548-I-3511 is well fitted by a 
soft component at low energy, and by a hard component 
above 2 keV. The X-ray flux increases of about 25 per cent 
in the two weeks between the XMM-Newton observations. 
The brightening is accompanied by a hardening of the 
spectrum, in particular above the energy break. This is a 
further indication that the emission in the high-energy part 
of the X-ray spectrum is dominated by the non-thermal jet 
emission, while a Seyfert component may be present in the 
low-energy part of the spectrum. 

• The broadband SED of both sources can be reproduced 
with synchrotron, SSG and external Compton scattering of 
the IR seed photons from the dust torus. The soft excess in 
J1548-I-3511 can be modelled as emission from a thermal 
corona. The disc luminosity of PKS 2004—447 turns out to 
be very weak. 

• The variability brightness temperature and the Doppler 
factor derived for both sources are similar to those found 
in 7 -ray blazars. These characteristics, together with the 


high radio-loudness are a good proxy for the presence of 
relativistic jet. However, they are not good tools in the 
selection of 7 -ray emitting NLSyl. 
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